Mon. Not. R. Astron. Soc. OOP. [HO ('2005') Printed 5 February 2008 (MN I*TeX style file v2.2) 



On the Rees-Sciama effect: maps and statistics 
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ABSTRACT 

Small maps of the Rees-Sciama (RS) effect are simulated by using an appropriate 
N-body code and a certain ray-tracing procedure. A method designed for the statis- 
tical analysis of cosmic microwave background (CMB) maps is applied to study the 
resulting simulations. These techniques, recently proposed -by our team- to consider 
lens deformations of the CMB, are adapted to deal with the RS effect. This effect and 
the deviations from Gaussianity associated to it seem to be too small to be detected 
in the near future. This conclusion follows from our estimation of both the RS angular 
power spectrum and the RS reduced n-direction correlation functions for n ^ 6. 

Key words: cosmic microwave background — cosmology: theory — large-scale struc- 
ture of the universe — methods: numerical — methods:N-body simulations. 



1 INTRODUCTION 

, From decoupling to present time, Cosmic Microwave Back- 
' ground (CMB) photons have been traveling -almost freely- 
' through a homogeneous and isotropic background universe 
with small perturbations. Taking into account the main re- 
sults from the analysis of the first year WMAP (Wilkin- 
, son Micr owave Anisot r opy P robe) data, which were pre- 
■ sented in lBennet et alj i2003fl . it is assumed: (i) a Gaussian 
distribution of adiabatic energy density fluctuations with 
Zel'dovich spectrum, (ii) a currently accelerated expansion 
with reduced Hubble constant h = 10~^Ho ~ 0.71 (Ho be- 
ing the Hubble constant in units of Km s~^AIpc~^), (iii) 
a cosmological constant with density parameter Ha = 0.73 
producing the current acceleration, and (iv) cold dark mat- 
ter (baryonic matter) with density parameter Qd ~ 0.23 
(f2f, = 0.04); hence, the total matter density parameter is 
Qrn = ilb + Hd = 0.27. In other words, a realization of 
the so-called conc ordance model is assumed and, then, the 
CMBFAST code dSeliak fc ZaldarriagallT99^ is used to get 
the linearly evolved power spectrum, P{k), of the scalar en- 
ergy density fluctuations, which is defined by the relation 
{6{k) S*{k')} = {2TTfP{k)S{k - k'). This spectrum is nor- 
malized with the conditio n gg = 0.93 to have the required 
density of galaxy clusters jEke. Cole fc FrenkllToQ^ . 

Along this paper, units are chosen in such a way that 
the speed of light is c = 1 and the gravitation constant is 
G = I/Stt. Whatever quantity "A" may be, A^, Ao, and 
stand for the A values on the last scattering surface, at 
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present time, and in the background, respectively. Symbols 
X, a, z, and rj, stand for the comoving position vector, the 
scale factor, the redshift, and the conformal time, respec- 
tively. Condition ao = 1 Mpc is assumed. 

Cosmological structures (scalar modes) produce a pe- 
culiar gravitational potential 0(5, 77) which interacts with 
the CMB. The partial time derivative of 0(a;, 77) produces 
temperature fluctuations. Quantity AT/T^ = {T -T^)/Tg 
measures fluctuations with respect to the CMB averaged 
temperature . This quantity is given by the formula: 



AT 

T ^ 



drj 



(1) 



where the integral is to be performed along the null geodesic 
of the perturbed universe associated to the observation di- 
rection n. As it is well known, a good approximation to 
AT/Tg can be obtained by calculating this integral along 
the background null geodesic corresponding to n, which 
obeys the equation 



— A(z)n* 
where 

X{z) = 



dh 



(2) 



(3) 



In the flat universe with cosmological constant we are con- 
sidering, linear, mildly non-linear, and strongly non-linear 
structures contribute to the integral in Eq. Q. Indeed, this 
integral allows us to estimate the anisotropy produced by the 
total time varying peculiar gravitational potential of the uni- 
verse. This anisotropy could be considered as an unique ef- 
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feet; nevertheless, the contributions -to the integral- due to 
linear and non-linear structure evolution are usually called 
the Integrated Sachs- Wolfe (ISW) and the Rees-Sciama (RS) 
effect, respectively. Some comments about this distinction 
seem to be worthwhile. During linear evolution, the pecu- 
liar potential can be written in the form (f)^D\{a) / a, where 
(f)^ is the initial potential (taken at decoupling time) and 
D\{a) is the linear growth function for the density contrast. 
At any time, the total peculiar potential, 0, can be splitted 
as follows: 



(4) 



where 



is the part of the potential properly due to non- 



linear evolution. After substitution in Eq. Q and integra- 
tion until present time, the first term of the r.h.s. of this 
equation gives the ISW effect, whereas the second one is the 
source of the RS effect. Equation I^J leads to a consistent 
definition of both effects (the addition is the total integrated 
gravitational efi'ect). Before the beginning of the nonlinear 
regime, function cft^^j^ takes on negligible second order val- 
ues and only the ISW anisotropy can be significant. During 
the nonlinear regime, the first term continues contributing 
to the ISW effect and the second term produces the RS 
anisotropy. Let us now study the ratio Di (a) /a to discuss in 
more detail ISW anisotropy production in the concordance 
model, in which, function D\ can be written as follows (in 
terms of z): 
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With the appropriate normalization factor, function Di is 
almost identical to the scale factor a for z > 3, whereas 
it is greater than this factor for 2: < 3. This behaviour is 
shown in Fig. Q, where we see that function Di/a is very 
close to (greater than) unity for z > 3 (z < 3). Two peri- 
ods can be then distinguished: in the first one {z > 3), the 
term (f)^D\{a)/ a is almost constant and no significant ISW 
anisotropy is generated. In the second period [z < 3), the 
potential cj>^D\(a) /a varies and the ISW effect is produced. 
These facts can be easily understood. For z > 3, vacuum 
energy is subdominant and the concordance model approx- 
imately evolves as a fiat universe without cosmological con- 
stant; in which, it is well known that no ISW is generated 
because the relation Di{a) ~ a is satisfied. All the ISW ef- 
fect is produced when vacuum energy plays a significant role 
in the evolution of the universe; namely, for z < 3; therefore, 
in the concordance model, the ISW and the RS effects are 
both produced at low redshifts. Only a small contribution 
to the RS effect is produced at z > 3 as a result of the 
first stages of non- linear evolution. The main part of the RS 
effect and all the ISW anisotropy are simulataneously gen- 
erated at redshifts z < 3. The ISW (RS) effect dominates 
for small (large) enough multipoles (see next sections). 

Since the location of the first peak of the CMB angular 
power spectrum indicates that the background universe is 
almost flat and, in such a case, the ISW effect only appears 
in the presence of vacuum energy (cosmological constant). 
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Figure 1. Ratio D\/a as a function of z in the concordance model 



any detection of the ISW effect could be considered as an 
indirect detection of vacuum energy. Since the largest ISW 
multipoles correspond to small I values with large cosmic 
variances, a direct detection based on the analysis of CMB 
maps becomes proble matic iHu fc Dod elson (200^ ; never- 
theless, ten years ago. ICrittenden fc Tur ok ( 199S) proposed 
an indirect detection method which is based on the estimate 
of cross-correlations between the CMB maps and tracers of 
the matter density distribution at z < 3 (risponsable for the 
ISW anisotropy). Many papers have been devoted to try 
this indirect detection jFosalba^Gaztanag^^^Castmid^ 



200; 
2004 



Fosalba fc GaztanagTil20M?lBm^^ 
Afshordi. Lho fc Straus;^ l'2004VNo ha et al.l J2004) : 



IVielva. Ma rtinez- Gonzalez fc Tucci (2006)); nevertheless, in 
the present paper, the ISW effect is only considered with the 
essential aim of estimating its contribution to our maps of 
the integrated gravitational effect, which are obtained from 
Eq. (0 and appropriate numerical techniques (see below). 
Since the ISW contribution is proved to be small enough (see 
next sections) , our maps can be considered as approximated 
simulations of the RS effect; by this reason, these maps are 
hereafter called RS maps and the spectra extracted from 
them are considered as an approximation to the RS spec- 
trum. 

This paper is organized as follows: Appropriate com- 
ments on the ISW effect are given in section 2. The nu- 
merical methods used to simulate maps of the RS effect are 
described in section 3. The RS angular power spectrum ob- 
tained from small maps is described in section 4, where the 
importance of the smallness of the maps (as a source of un- 
certainty) is discussed in detail. Next section deals with the 
estimation of both reduced m-correlations and deviations 
from Gaussianity and, finally, section 6 contains a summary 
of conclusions and a general discussion. 



2 ON THE ANGULAR POWER SPECTRUM 
OF THE ISW EFFECT 

In a flat background with cosmological constant, the angular 
power spectrum of the ISW effect depends on various param- 
eters. The values of h, S7a, and as are particularly impor- 
tant. This spectru m has been previ o usly calculated in var - 
lous cases, see, e.g. iHu fc Dodeisoiil (l2002h : ICooravl J2002l) . 
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Figure 2. ISW power spectrum, A^, in fiK, as a function of 
\og(£). 
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Figure 3. rms relative mass fluctuations inside a sphere whose 
present radius in Megaparsecs is R. 



Here, it is calculated for the cosmological model and param- 
eters fixed in tfie introduction. T lie Ci quantities of tli e fSW 
effect can be written as follows jFuUana fc Sae 



Ce = N 



where 



P{k) ,2 



M[\k]—[{l + z)Di{z)]d\ 



ji being the spherical Bessel function of order I, and 



N = 



Dm 



(7) 



(8) 



(9) 



where function D\(z) is given by Eqs. @ and ©. The in- 
tegrals in Eqs. JTJ and Q have been numerically calculated 
to get the required spectrum. Results are presented in Fig. 
©, where quantity A^, = + l)C«/27r]^/^ in (used 
along the paper to describe the angular power spectrum) is 
displayed in the I interval [2,700]. 

Our estimation, as we ll as pre vious ones (see e.g. 
iHu &, Dodeisoiill|2002l1 : ICoor av ( 20021 : |Padmanabhan etall 
l|2005|v r. leads to two main conclusions: (1) the most impor- 
tant part of the ISW effect is produced by very large spatial 
scales (small k values) contributing to small multipoles, and 
(2) the contribution to this effect strongly decreases as the 
spatial scale does. For the parameters we have assumed, the 
A.J, value corresponding to ^ = 700 is 0.037 and, more- 
over, quantity A^, rapidly decreases as I increases reaching 
a value close to 0.0045 for I ~ 2000. These values are 
to be compared with those extracted from our numerically 
simulated maps. As a technical requirement associated to 
ray-tracing (see below), our maps are produced by the total 
time varying peculiar potential due to all the scales smaller 
than 60 Mpc, which are evolved by N-body simulations from 
redshift 5.2 to present time; hence, all the ISW effect pro- 
duced (at 2: < 3, see Fig. Q) by these scales is included 
in the resulting maps. The effect produced by scales greater 
than 60 Mpc must be independently calculated. Is it a pure 
ISW effect? A few considerations about linearity are useful 
to answer this question. 

In order to discuss about linearity in detail, the root 
mean square (rms) of the relative mass fluctuations inside a 



sphere having a present radius of R Mpc (randomly placed 
in our flat universe) is estimated. Condition ao = 1 implies 
that number R is also the comoving radius of the sphere. 
This rms value is g{R) = (AM/M)™, = {{SM/MfY^"^, 
with 



{{5M/MY) 



1 

2^ 



k^P{k)W^{kR)dk 



(10) 



where W{kR) is the window function of the R-sphere 
[W{kR) = -\{smy — ycosy) with y = kR; see IPeeble j 
(1980)]. Numerical calculations based on these formulae lead 
to the a{R) values presented in Fig. |3] From this figure it 
follows that, for the chosen normalization of P{k) and the 
comoving radius i? = 15 Mpc, the (715 value is close to 0.7; 
hence, it can be assumed (standard point of view) that re- 
gions having a comoving diameter greater than 30 Mpc can 
be treated as linearly evolving zones. The scale L = 60 Mpc, 
used along this paper, evolves in the linear regime with 
(t(30) ~ 0.4; hence, scales greater than 60 Mpc would pro- 
duce a pure ISW effect to be separately estimated with ap- 
propriate Unear techniques. The ISW effect produced by all 
the scales smaller than 60 Mpc is included in our maps and, 
obviously, it is smaller than the total ISW effect (see Fig. |5| 
and previous comments). 



3 NUMERICAL METHODS 

RS is a pure gravitational effect and, consequently, only the 
dominant dark matter component is considered, whereas the 
sub-dominant baryonic component is neglected. On account 
of these facts, the formation and evolution of non-linear 
cosmological structures can be described by using N-body 
simulations with appropriate boxes and resolutions. In the 
PM simulations we use he re jHocknev fc Eastwoodl [l983 : 
iQuilis. Iba nez fc Sae the peculiar potential satisfies 

the equation ( Martinez- Gonzalez. Sanz fc Silk . 1990.') ; 



A<^ = -a (pm — Pms) 



(11) 



The universe is assumed to be covered by simulation boxes 
and, consequently, it is periodic. A s it was pointed out b y 
ICerda-Duran. Ouihs fc Sae j ^|2004^ and lAnton et al.N2005l) . 
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periodicity effects magnify lens deformations and, in gen- 
eral, gravitational anisotropies as the RS effect. Various 
techniques (ray-tracing methods) h ave been used to avoid 
this m agnification, for example, iTuluie. Laeuna fc Anniiio^ 
il996i) averaged the temperature contrasts of many rays, 
but then, they found a map with too large pixels, from 
which, the multipoles corresponding to ^ > 700 could not 
be obtained. Another method based on multiple plane pro- 
jections was proposed and ap plied to study week l ensing 
by cosmological structures (see ljain. Seliak fc White 
for a detailed description). Afterwards. IWhit^^^ 
designed the "tiling" ray tracing procedure, in which, inde- 
pendent simulations with appropriate boxes and resolutions 
tile the photon trajectories. Finally, another method based 
on the existence of preferred directions and on the use of 
an appropriate cutoff was recently proposed and applied 
fCerda-Duran. Quilis fc Sacz (2004); Anton ct al. (200^). 
In the required cutoff, the Fourier modes corresponding to 
spatial scales larger than a given one (Lmax) are eliminated 
from the peculiar gravitational potential. The Lmax value 
is chosen in such a way that: (a) the cutoff eliminates all 
the spatial scales which are too large to be well described in 
the simulation box, and (b) after the cutoff is performed, 
CMB photons cross neighbouring boxes through statisti- 
cally independent regions. This cutoff is only performed in 
the output peculiar gravitational potential with the essen- 
tial aim of calculating the integral in Eq. Q; however, all 
the spatial scales allowed by the box size are taken into ac- 
count in the N-body simulation. Our ray-tracing method is 
a very good choice pres enting some advanta ges with respect 
to other techniques (see lAnton et al.l ('2005)), e.g., there are 
no discontinuities at the crossing points between successive 
boxes, the computational cost is moderated, and so on. 

In order to study the RS effect, we devel op the same 
type o f N-body simulations which were used, in I Anton et alJ 
(|2005|), to deal with lensed maps of the CMB. In this way, 
various maps of the relative temperature variations due to 
the RS effect are created and analyzed, and the dependence 
of the results on various parameters involved in the simula- 
tions and also in their statistical analysis is then discussed. 
Most maps are regularly pixelised, but some maps based on 
HEALPIx (Hiera rchical Equal Area Isolatitude Pi xelisation 
of the Sphere, see iGorski. Hivon fc Wandell (ll999M are also 
considered. Taking into account that t he box size used in 
the mentioned N-body simulations fsee lAnton et alJ ll2005ll 
for details) is L = 256 Mpc, and also that the main part of 
the RS effect is assumed to be produced at redshift 2 < 5.2 
(see below), the size of the resulting RS maps appears to 
be close to 2°. Simulation boxes cover a periodic universe 
and, for the direction 9 = 77.2° and = 12.6° (where 6 
and (j) are spherical coordinates defined with respect to the 
box edges), the CMB photons cross around 30 boxes be- 
fore passing through the region where they were initially 
located (in the first box) . It can be verified that these boxes 
cover the trajectory of the CMB photons from the redshift 
Zin ~ 5.2 to present time, and also that these photons cross 
successive boxes through independent regions separated by 
a distance close to 52 Mpc. The cutoff is performed at the 
scale kmin = 2-K/Lmax with Lmax = 60 Mpc (in agreement 
with previous comments about the scales involved in our 
simulations) . 
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Figure 4. Simulated 1.83° X 1.83° map of the RS effect obtained 
from a HR N-body simulation 



Along the paper we develop three types of N-body sim- 
ulations in boxes of L Mpc (integer and even L): Low Res- 
olution (LR) simulations involving L/2 cells per edge and 
(L/2)'' particles (2 Mpc cell size and ~ 380 time steps). 
Intermediate Resolution (IR) simulations using L cells per 
edge and particles (1 Mpc cell size and ~ 410 time steps), 
and High Resolution (HR) simulations with 2L cells per edge 
and (2L)^ particles (0.5 Mpc cell size and ~ 500 time steps). 
Figure I^J shows a RS map which has been obtained with 
our ray-tracing method and a HR simulation of structure. 
Its look is promising, but the correlations appearing in this 
type of maps must be estimated to get objective conclusions. 
It is done in next sections. 

For any observation direction n, function d(j>/dri must 
be calculated in a set of points to estimate the integral in Eq. 

We have chosen equidistant points (on the correspond- 
ing background geodesic) separated by a comoving distance 
equal to the cell size. Each of these points has a comov- 
ing position vector Xp , and the photon passes by this posi- 
tion at time rjp . Point P is placed inside a certain P-cell of 
the computational grid and time rjp belongs to some inter- 
val [rii,r]i+i], where rji and r/i+i are two successive times of 
the N-body simulation, at which, the peculiar potential (j) is 
found at any grill node. The potential at point P is calcu- 
lated at times rji-i, rji, rji+i, and ?7i+2. The calculation at 
any of these times is performed by means of a suitable 3D 
interpolation, which uses the potential in the vertices of the 
P-cell and in those of the neighbouring ones. From the re- 
sulting potential at the above four times, the partial deriva- 
tive d(j)/drj is finally calculated with an appropriated numer- 
ical method. For each type of N-body simulations (LR, IR 
or HR), various methods for spatial interpolation and also 
for the numerical calculation of time derivatives have been 
tried. The final results do not significantly depend on the 
chosen methods because these results depend much more on 
the simulation resolution (see comments on Fig. ^ given 
below). Two types of spatial interpolations have been con- 
sidered: (i) an standard 3D linear interpolation using the 
potential in the eight nodes of the P cell, and (ii) a 3D in- 
terpolation based on the use of splines along each direction; 
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Figure 5. RS power spectrum, A^, , in fiK, as a function of log{£). 
Solid, dotted and dashed lines are the spectra produced by spatial 
scales smaller than 60 Mpc, 30 Mpc, and 15 Mpc, respectively. 
IR simulations have been used. 
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Figure 6. RS power spectrum, A , in fiK, as a function of 
log{£). Solid, dashed and dotted lines are the spectra obtained 
from the LR, IR, and HR simulations defined in the text, respec- 
tively. Scales are smaller than 60 Mpc in all cases. 



in this case, four vertices located in the P-ceU and in their 
neighbouring ones are considered for interpolations in the 
directions of the box edges; results from both interpolations 
are very similar; relative differences between the Ce coeffi- 
cients extracted from the maps obtained with the two inter- 
polations are smaller than 0.02 for any £. Also time deriva- 
tives based on different numbers of points, rji, have been 
performed, but the results are almost identical in all cases 
because </> is a very smooth function of rj in any node of the 
simulation grid (^see lAliaea. Quilis fc SiAe2 f2002")). In short, 
a 3D linear interpolation in space and a standard derivative 
based on two times r/i and rji+i suffice to get good maps 
of the RS effect. Various methods have been also used to 
perform the integral in Eq. Q: trapezoidal, Simpson, and 
so on. Results are almost identical in all cases. All these 
comments show that our ray-tracing procedure is a robust 
numerical method, which leads to RS maps having almost 
the same properties for a wide set of numerical techniques 
(interpolations, derivatives and integrations). 



4 ON THE ANGULAR POWER SPECTRUM 
OF THE RS EFFECT 

Since we are concerned with small squared regularly pix- 
elised maps, the power spectrum estimator described in 
lArnau. Aliaga fc^ acz (2002) applies. Although this estima- 
tor is used as the basic one along the paper, some results 
are compared with those of the code ANAFAST included in 
the HEALPIx package. 

The basic estimator is first used to analyze 1.83° x 1.83° 
RS maps obtained from IR simulations with L — 256 Mpc. 
Figure |S] shows the resulting angular power spectra (A^) for 
three appropriate choices of the cutoff: 60, 30, and 15 Mpc. 
Differences between the values corresponding to the 
solid (dotted) and dotted (dashed) lines are produced by 
scales between 60 and 30 Mpc (30 and 15 Mpc). The rela- 
tion Aj, > 0.037 is satisfied and, according to previous 
comments (see section 2), the ISW effect cannot produce 
these spectra; therefore, scales between 60 Mpc and 15 Mpc 
seem to produce a certain Rees-Sciama effect, which is only 



possible if they become non-linear enough in our simula- 
tion boxes. The question is: are the scales between 15 and 
60 Mpc well described in our simulations? In order to an- 
swer, it should be taken into account that the evolution of 
the scales under consideration depend on N-body simula- 
tion characteristics as resolution and box size; indeed, the 
greater the N-body resolution, the better (more realistic) 
the description of all the spatial scales and, in particular, 
the description of the scales between 15 and 60 Mpc. We 
can say then that part of the non-linear behaviour of these 
scales is due to low resolution (see next paragraph for more 
discussion) . 

For a cutoff at 60 Mpc and a fixed box size L — 
256 Mpc, the angular power spectrum of the RS maps has 
been calculated for HR, IR, and LR simulations; namely, for 
different space and mass resolutions. The resulting spectra 
are displayed in Fig. |^ where we see that the greater the 
resolution the smaller the amplitude of the RS effect. Tak- 
ing into account Fig. |K| this behaviour of the amplitude is 
not surprising because, as resolution increases, the evolution 
of the scales between 60 and 15 Mpc is better described and 
masses concentrate more and more on smaller and smaller 
scales. Taking into account that the separation between the 
solid and dashed lines is greater than that of the dashed and 
dotted lines, it seems that, although the spectra of Fig.|S|de- 
pend on resolution, they tend to a certain spectrum located 
below the dotted line and close to it; hence, this line gives 
an upper limit to the true spectrum. The difference between 
the true spectrum and the upper limit (dotted line) is ex- 
pected to be small and, consequently, the dotted line can be 
considered as our best estimation of the RS angular power 
spectrum. In spite of this fact, there are systematic errors in 
this spectrum associated to the smallness of our simulated 
maps. These errors are analyzed below 

The spectra of Fig. |S| have been obtained for an ini- 
tial redshift Zi„ = 5.2, but other initial redshifts have been 
tried. Fig.|7|shows the spectra obtained from initial redshifts 
of 5.2, 10.2 and 15.2. In the three cases, IR simulations have 
been used and the size of the simulated maps is 1.83° x 1.83°. 
In order to maintain the map area for the three initial red- 
shifts, the size of the simulation box (and accordingly the 
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Figure 7. RS power spectrum, A^, , in ^K, as a function of log(£). 
Solid, dashed, and dotted lines correspond to the spectra obtained 
from the initial redshifts 5.2, 10.2, and 15.2, respectively. IR sim- 
ulations and scales smaller than 60 Mpc have been considered in 
the three cases. 



resolution) has been varied. The assumed box size (N-body 
resolution) slightly increases (decreases) as the initial red- 
shift grows. Figure |7| shows oscillations in whose am- 
plitude increases as Zin does. In spite of these unexplained 
oscillations, the resulting spectra are rather similar in the 
three cases and, as it follows from this figure, the choice of 
the initial redshift Zin ~ 5.2 is the simplest one leading to a 
good estimate of without spurious oscillations. 

The spectra obtained from our small non-Gaussian RS 
maps have systematic errors which are now approximately 
estimated. In order to perform this estimation, we begin with 
the analysis of Gaussian maps of the same size and, then, 
results are considered as an approximated estimate of the 
errors associated to the non-Gaussian RS simulations. This 
procedure seems to be reasonable taking into account that 
deviations from Gaussianity of the RS maps do not seem to 
be dramatic (see our discussion in section 5). As it is well 
known, the relative errors of the Ce quantities extracted 
from a Gaussian map is ACi/d = [2/{2£ + l)Y^^f~kJ^, 
where fsky is the ratio between the map area and th at of 
the fuU sky llScott. Sredniki fc Whitet l994: Knt^HH- Ac- 
cording to this formula, the relative errors corresponding to 
3.66° X 3.66° Gaussian maps range between ~ 2.5 and ^ 0.5 
as £ takes on the values of the interval (500, 10000). For 
1.83° X 1.83° (7.2° X 7.2° ) maps these errors increase (de- 
crease) by a factor 2. These relative errors are so large that 
one could think that Gaussian maps of dominant anisotropy 
(and also our simulated maps) are not useful at all; how- 
ever, when C'l quantities are extracted from one of the men- 
tioned Gaussian maps, the errors lead to a oscillatory an- 
gular spectrum which changes from £ to £ + 1; in other 
words, the resulting spectrum has high frequency oscilla- 
tions which can be eliminated by means of suitable averages. 
These averages are performed during (after) the Ce estima- 
tions when the method proposed in lArnau. Aliaea fc Saej 
120021 (ANAFAST) is used. The final averaged spectra are 
meaningful in a certain interval of £ values. In order to de- 
fine this interval, 1.83° x 1.83° maps of the CMB domi- 
nant anisotropy have been simulated by using: the angu- 
lar power spectrum given by CMBFAST (in the model as- 
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Figure 8. CMB power spectrum, A^, in fiK, as a function of 
\og{£) in all panels. All the solid lines exhibit the CMB spectrum 
given by CMBFAST which is used as an input to simulate CMB 
maps. Top, middle, and bottom panels correspond to 1.83° X 
1.83°, 3.66° X 3.66°, and 7.2° X 7.2° maps, respectively. In the 
three panels, dotted and dashed lines are the spectra extracted 
from the simulated maps by using two distinct estimators (see 
text). 

sumed here) and a map making method based on the fast 
Fourier transform Sacz. Holtmann & Snioot (199^. The re- 
sulting maps are analyzed with two different power spec- 
trum estimators. In the top panel of Fig. |Sl the solid line 
is the CMBFAST spectrum used as input in the simulated 
maps, the dashed line is the spectrum extracted from the 
maps by using the p ower spectrum estimator described in 
lArnau. Aliaga fc Sa cz ( 2002) and, finally, the dotted line has 
been obtained using ANAFAST (see above) to analyze an 
equatorial healpix map obtained from our uniformly pix- 
elised one. Both estimators give very similar spectra which 
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Figure 9. RS power spectrum, A^,, in fiK, as a function of 
log{£). Dotted-dashed and dashed (solid and dotted) lines have 
been obtained with LR (IR) simulations for L = 256 Mpc and 
L = 512 Mpc, respectively. The angular sizes of the maps are 
given inside the panel. Scales are smaller than 60 Mpc in all 
cases. 



can be considered as a rather good estimate of the true spec- 
trum for 500 ^ ^ 2000. Maps are too small to get the 
spectrum for £ < 500 and, furthermore, for £ > 2000, the 
signal is too weak and, then, numerical and statistical noise 
become comparable to it. The middle and bottom panels 
of Fig. |S] have the same structure as the top panel; how- 
ever, in the middle (bottom) panel the angular size of the 
maps is 3.66° x 3.66° (7.2° x 7.2°). Comparing the panels 
of Fig. |H1 one concludes that the more extended the maps, 
the better the resulting spectra and the wider the £ interval 
where the spectra are good enough. In the case of the RS 
effect studied with HR simulations, the size of our simula- 
tion boxes (L = 256 Mpc) and, consequently, the angular 
size of the resulting maps (1.83° x 1.83°) cannot be freely 
increased because any significant increment would lead to a 
problematic decrease of the N-body simulation resolution, 
nevertheless, if the above results corresponding to Gaussian 
maps are applied -eis an approximation- in the presence of 
deviations from Gaussianity, our 1.83° x 1.83° maps should 
suffice to get a rough, but useful, estimate of the RS spec- 
trum in a wide £ interval. This conclusion is confirmed by a 
direct study of RS maps presented in the next paragraph. 

In order to have direct information about the qual- 
ity of the spectra associated to 1.83° x 1.83° RS maps, 
two tests have been performed. In the first one, maps of 
this angular size have been obtained from LR simulations 
with L = 256 Mpc and, then, the angular power spec- 
trum extracted from these maps (which is presented in the 
dashed-dotted line of Fig. ^ is compared with that ob- 
tained from LR simulations with L — 512 AIpc which lead 
to 3.66° X 3.66° RS maps (dashed line of the same Fig- 
ure). Since both simulations have the same resolution, dif- 
ferences in the resulting spectra are due to the difference 
between the angular size of the maps. The comparison of 
the dashed-dotted and dashed lines of Fig. |n| shows that, in 
the 1.83° X 1.83° maps, the spectrum is slightly underesti- 
mated for 500 <£< 800 and overestimated for £ > 800. The 
second test is identical to the first one except for the reso- 
lution of the involved simulations. The 1.83° x 1.83° maps 



are obtained with IR simulations and L — 256 Mpc, and 
the maps whose angular size is 3.66° x 3.66° are created 
with IR simulations in boxes with L — 512 Mpc; the cor- 
responding spectra are displayed by the solid and dotted 
lines of Fig. |5| respectively. The comparison of these spectra 
indicates that, for the smallest maps, the resulting spec- 
trum is underestimated (overestimated) for 500 < £ < 700 
(£ > 700). Furthermore the maxima of the dashed-dotted 
and solid lines appear as a result of the overestimations and 
they disappear in the case of our most extended maps (see 
the dotted and dashed lines) . Unfortunately, HR simulations 
cannot be used to perform a third test similar to the previous 
ones (too great computational cost of HR simulations with 
L = 512 AIpc), nevertheless, previous discussion strongly 
suggests that, also in the case of 1.83° x 1.83° RS maps gen- 
erated with HR simulations, the spectrum (dotted line of 
Fig. |SJ should be underestimated and overestimated as in 
the first and second tests (see the weak maximum observed 
in that dotted line); therefore, for £ > 800, the A^, quantities 
corresponding to 3.66° x 3.66° HR maps should be smaller 
than those extracted from 1.83° x 1.83° HR maps (dotted 
line of Fig. |SJ. In other words, the smallness of the maps 
leads to an overestimation of the angular power spectrum 
for large £ values. 

Previous discussion about the uncertainties associ- 
ated to both N-body resolution (Fig. and map size (Fig. 
|UJ can be summarized as follows: (1) for £ values between 700 
and 800, the spectrum A^ slightly depends on the map size 
(see Fig. 1^; however, its dependence on N-body resolution 
(see Fig.lSJ is important; hence, the A^ value corresponding 
to ^ ~ 700 essentially depends on resolution. From the visual 
analysis of Fig.|S|one easily concludes that this value should 
be a number between 0.4 jiK and 0.5 nK, (2) for £ > 800, 
the overestimations of A^ associated to N-body resolution 
and map size are comparable and they add to generate a 
total excess of power which might be close to 0.2 jiK. This 
number follows from a qualitative study of Figs. |S| and 1^1 
(3) on account of points (1) and (2), the dotted line of Fig. 
|S] appears to be an upper limit to the RS spectrum, but it 
can be considered, in practice, as an useful rough estimation 
of the true spectrum. The most important point is that the 
resulting spectrum has been obtained from fully nonlinear 
dark matter evolution (without approximating techniques o 
models), and (4) more resolution and more extended maps 
are necessary to get a more accurate RS spectrum. 



5 STATISTICAL ANALYSIS 

Various reduced correlation functions are estimated and 
then deviations from Gaussianity are analyzed. The corre- 
lations associated to n directions (or n points on the last 
scattering surface) are estimated for n ^ 6. Given a map of 
the variable and m directions, the reduced angular corre- 
lations are: 



Cm = {C(ni)C(^2)---C(n™)) 



(12) 



where the average is over many statistical realizations. 

Correlations depend on the relative positions of the cho- 
sen directions; hence, for CMB and RS maps, they depend 
on the figures that these directions draw on the last scatter- 
ing surface. In this paper we have used the sets of directions 



8 N. Puchades et al. 



Figure 10. m-direction configurations for statistical analysis 

displayed in Fig. IIUI where the basic correlation angle, a, is 
that subtended by the two directions of the case w = 2 . In 
Gaussian statistics it is well known that llPeeblej|l98Cl) : (i) 
all the Cm correlations vanish for odd m and, (ii) for even m, 
all the correlations can be written in terms of C2. Hereafter, 
these even Gaussian correlations are denoted Cgm (the suf- 
fix g stands for Gaussian) . For the configurations of Fig. 1101 
one easily obtains the relation: 

CgA{a) = 2C|(q) + Cl{V2a) , (13) 

leading to Cg4(0) = 3(71(0) for a = 0, and the equation: 

Cg(i[a) = 3C|(a) + 2C2(a)C|(%/2a) 

+ 4C2(a)C2(2Q)C2(V2a) 

+ 2C2(V5q)C|(\/2q) + C2(a)C|(\/5a) 

+ 2C2(%/5q)C|(q) + 2C2(a)C|(2«) , (14) 

which takes on the form Cg&{Q) = 15C|(0) for a = 0. Our RS 
simulated maps are not Gaussian and we are interested in 
deviations from Gaussianity. Even correlations can be used 
to estimate these deviations as follows: (a) extract the C2, 
C4, and Ce correlations from the simulated maps, (b) use the 
extracted C2 correlation to calculate Cg4, and Cge functions 
assuming Gaussianity, namely, using Eqs. (1131 - 1141 1 and, (c) 
compare C4 and Ce correlations with functions Cgi and Cg6, 
respectively. Differences between C4 and Cg4 as well as be- 
tween Ce and Cg6 measure deviations from Gaussianity. Fur- 
thermore, odd correlations vanish in the Gaussian case and, 
consequently, if at least one of the correlations C3 and C5 
(for the configurations of Fig. IIOII appears to be different 
from zero, we can conclude that the analyzed RS maps are 
not Gaussian. 

In order to get the Cm correlation from a set of maps, 
the m-configuration of Fig. 1101 is randomly placed on the 
maps and, then, the temperature contrasts in the nodes of 
the configuration are calculated by appropriate interpola- 
tions in our maps (making use of their pixelisation) ; finally, 
the average of Eq. p2|l is calculated considering all the maps 
and all the locations of the m-configuration. As it occurs 
in the case of the Ce quantities, the fact that our simu- 
lated maps are small implies that the Cm correlations ex- 
tracted from them are not the true ones. In spite of this 
fact, it is now proved that the methods proposed above to 
detect deviations from Gaussianity -based on the use of even 
correlations- work. In order to do that, these methods are 
first applied to CMB maps -Gaussian by construction- hav- 
ing either the same sizes as our simulated RS maps or larger 
ones; in this way, small Gaussian maps are characterized 
and, then, deviations of the RS maps with respect to this 
characterization point out a certain level of non- Gaussianity. 

We begin with 1.83° x 1.83° Gaussian CMB maps. Some 
correlations extracted from these maps are displayed in the 
panels of Fig. 1111 The dashed (dashed-dotted) line of the top 
panel exhibits the C4 (Cg4) correlations (see above), whereas 
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Figure 11. Even correlation functions obtained from 1.83° X 
1.83° maps in terms of the correlation angle a. C4 (Ce) correla- 
tions are multiplied by an appropriate power of ten and given in 
fiK'^ (fiK^). Top (bottom) panel corresponds to C4 ( Ce ) func- 
tions. Dashed-dotted (dashed) line exhibits either Cg4 or Cg6 ( 
C4 or Ce ) correlations extracted from CMB maps. Solid (dotted) 
line shows either Cg4 or Cg6 ( C4 or Ce ) correlations found from 
RS maps. Each type of line is defined inside the panel giving: (i) 
the power of ten exponent, (ii) the type of map (CMB or RS), 
and (iii) the type of correlation (S, for C4 and Ce, and G, for Cg4 
and Cge) 



the same lines of the bottom panel display Ce and Cge cor- 
relations. Since the 1.83° x 1.83° maps are too small, cor- 
relations extracted from them are not the true correlations 
(which could be obtained from a large set of big maps) and, 
consequently, in spite of Gaussianity, functions C4 and Cg4 
slightly separate (compare dashed and dashed-dotted lines 
in the top panel of Fig. lllH . The same occurs for the correla- 
tions Ce and Cge (compare dashed and dashed-dotted lines 
in the bottom panel of Fig. llH . In the case of 3.66° x 3.66° 
maps, the correlations are presented in Fig. 1121 using the 
same notation as in Fig. 1111 From Fig. 1121 one easily con- 
cludes that dashed and dashed dotted lines are much closer 
than in the corresponding panels of Fig. 1111 It has been ver- 
ified that for 7.2° x 7.2° maps, the corresponding lines are 
closer than in Fig. 1121 This behaviour characterize Gaussian 
maps, in which the C4 and Ce correlations approach the Cg4 
and Cge functions, respectively, as the map size increases 
and the extracted correlations approach more and more the 
true ones. It is worthwhile to notice that, although the cor- 
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Figure 12. Same as in Fig.EJfor 3.66° X 3.66° maps 



relations C2, C4 and Ce obtained from small maps are not 
the true correlations, they allow us to verify the Gaussian 
behaviour described above, which can be applied to detect 
deviations from Gaussianity. 

Even correlations extracted from 1.83° x 1.83° and 
3.66° X 3.66° RS maps are now analyzed. In order to do 
that, the study described in the last paragraph (for CMB 
maps) is repeated for RS maps. Results are also presented 
in Figs. llll and ll2l in which solid and dotted lines are used to 
represent the same correlations given by the dashed-dotted 
and dashed lines of the same Figures (see above) . From the 
comparison of the solid and dotted lines one easily concludes 
that: (i) for 1.83° x 1.83° maps fFig. Illfl. these lines sepa- 
rate for 0.07° < 9 < 0.15°. Since this separation is compa- 
rable to that corresponding to Gaussian CMB small maps 
(dashed-dotted and dashed lines), it does not point out any 
deviation from Gaussianity, and (ii) for 3.66° x 3.66° maps 
fFig. Wli) . the separation between solid and dotted lines is 
much greater than that of the CMB Gaussian maps with 
the same size (dashed dotted and dashed lines of Fig. 112^ . It 
occurs for all the angular scales we have represented, which 
include the scales corresponding to 500 < i < 2000, for 
which ou best estimation of the RS spectrum (dotted line 
of Fig. |SJ has been obtained. This result (great separation) 
proves that the RS maps are not Gaussian because, in the 
case of Gaussianity, the separation of the solid and dotted 
curves corresponding to 3.66° x 3.66° maps must be smaller 



than the separation corresponding to 1.83° x 1.83° maps 
(Gaussian behaviour defined above). 

Odd correlations have been also considered. Maps of 
1.83° X 1.83° and 3.66° x 3.66° squared degrees have been 
used to study these correlations. In the case of CMB Gaus- 
sian simulations, the general method (see above) for esti- 
mating functions C3 and C5 has been applied using n small 
CMB maps. Different n values have been considered. Results 
do not converge as n increases because the true odd corre- 
lations vanish and, consequently, the functions C3 and C5 
extracted from n simulated Gaussian CMB maps are pure 
errors for any n. This complete absence of convergence char- 
acterizes small Gaussian maps for any size. In the case of RS 
small maps with any of the above sizes, it has been verified 
that the Cs and C5 correlations extracted from a number m 
of these maps do not converge as m-increases (as it occurs 
for Gaussian maps); this means that we have not been able 
to detect the small odd correlations associated to the RS 
effect (see next section for more discussion). 



6 RESULTS AND DISCUSSION 

In this paper, an improved version of the numerical 
methods used in Aliaga, Quilis & Sacz (2002) is applied 
to study the RS effect. Our main improvements: pre- 
ferred direction and appropriate cutoff were described in 
ICerdir^DuraJij(^ ^ Saed (|20o3) and successfully applied 
ii n^it6i^^alTll200^) ~These new elements make negligible 
the magnification of the RS and lens effects due to the spa- 
tial periodicity of the simulated universe, which is covered by 
identical simulation boxes. The present size of these boxes 
must range between ~ 256 Mpc and ~ 512 Mpc. Only in 
this way, good N-body simulations with large enough reso- 
lutions can be obtained using our PM code; as a result of 
this size restriction, we are constrained to work with small 
maps of the RS effect having sizes between 1.83° x 1.83° 
and 3.66° x 3.66° (see section 3). Since the correlations ob- 
tained from these small maps are expected to have impor- 
tant uncertainties, an exhaustive study about the accuracy 
and significance of these correlations is developed. Further- 
more, methods specially designed for detecting deviations 
with respect to Gaussianity (using the correlations of small 
maps) are described and applied. 

The analysis of small CMB Gaussian maps of dominant 
anisotropy (CMBFAST angular power spectrum, see section 
4) has been basic in order to study both power spectrum esti- 
mations and non-Gaussianity detection. The following stud- 
ies have been developed: (1) the Ce quantities extracted from 
Gaussian CMB maps (for different sizes) have been com- 
pared with those used in their simulation; in this way, the 
interval {£1/2) where the recovered angular power spectrum 
is good enough (from a qualitative point of view) has been 
estimated (see Fig. |S]of section 4), the existence of this in- 
terval has been interpreted taking into account the concept 
of sample variance and, (2) the method described in section 
5 to look for deviations with respect to Gaussianity -from 
even correlations- has been applied to analyze CMB maps 
with distinct sizes and, thus, Gaussian maps have been char- 
acterized, which has been basic for further researches about 
deviations from Gaussianity. The interval (^1,^2) depends on 
the map size (see section 4) and the Ce quantities have been 



10 N. Puchades et al. 



calculated for various of these sizes. The larger the size of 
the RS map, the wider this interval. More extended maps 
are necessary to get the angular power spectrum for £ < £i, 
and more extended maps and greater resolutions are needed 
to find this spectrum in the case £ > £2', namely, for very 
small angular scales. There is the opinion that the RS effect 
could produce some feature in the angular power spectrum 
at small enough angular scales, by this reason, our main 
project in this field is the use of P3M and AP3M N-body 
codes to get resolutions as high as possible in large enough 
boxes. In spite of this plan, we are skeptical about the exis- 
tence of the mentioned feature because the RS effect is pro- 
duced by the time derivative of the peculiar potential, which 
seems to be rather well estimated with the effective resolu- 
tion of 1 Mpc corresponding to the dotted line of Fig. 15] 
in which, there is no any trace of incipient features at large 
£ values; nevertheless, it must be recognized that, for these 
large values, the uncertainties in our spectrum are maximum 
(superimposition of the overestimations due to resolution 
and map smallness) and, consequently, a moderated excess 
of power due to fully non-linear evolution could be hidden 
by our current uncertainties. Our method (based on N-body 
simulations) is fully nonlinear and, moreover, it takes into 
account both the effect due to the gravitational collapse of 
cosmological structures and that produced by the motion of 
these structures in the box. These are the most promising as- 
pects of the proposed techniques whose future applications 
-based on better N-body simulations and computers- can 
improve on both the estimation of the RS spectrum (point- 
ing out posible new features) and the detection of deviations 
from Gaussianity. In other words, although improvements on 
the applications of the ray tracing procedure used here are 
possible, current applications lead to too small RS maps and 
moderate N-body resolutions and, consequently, only an up- 
per limit to the RS spectrum (rough estimate of it) has been 
obtained in this paper. Fortunately, we are going to see that 
this estimate suffices to discuss detection of the RS effect 
with WMAP and PLANCK sateUites (see below) and, fur- 
thermore, our detection of deviations from Gaussianity in 
RS maps seems to be doubtless. 

A few comments about previous results are 
now necessary for comparisons. The angular power 
spectrum of the Rees-S ci ama effect was studied by 
iTuluie. Laeuna fc Anniiio^ Jl99d) and ISeliakI Jl99et) . In 
both papers, the effect was estimated for open and flat 
cold dark matter universes without cosmological constant. 
In the first of these papers, the authors used N-body 
simulations and a certain ray-tracing procedure to get 
the Ce coefficients for £ < 700. Among their results we 
emphasize the following ones: (a) in the case of a flat 
universe (with h=0.75 and as = 1), and for £ close to 700, 
quantity takes on values around 0.27 fiK, and (/3) for 
£ ~ 100 this quantity reac hes a maxim um value close to 
0.85. In the second paper, ISeliakI lll99d) did not use any 
ray tracing procedure. He used N-body simulations (second 
order perturbation theory) to evolve the density contrast 
of strongly (mildly) non-linear structures and, then, from 
the resulting contrast, the author estimated the power 
spectrum of the quantity d(j){x,t) / dt. Finally, from this 
spectrum, the Ci quantities were calculated. Four cases 
corresponding to different values of the parameters erg and 
Q,mh were studied using N-body simulations. Some results 



corresponding to these four cases are now summarized: (i) 
the temperature contrast (AT/T) takes on values between 
10"'' and lO"'^, (ii) quantity is maximum between 
£ = 100 and £ = 300, (iii) the maximum values of A^ range 
from 0.38 ^j.K to 2.25 fj.K and, (iv) the values corresponding 
to ^ ~ 700 ranges from 0.27 fiK to 1.2 fj,K (the first of 
these values corresponds to a flat universe with h=0.5 and 
(Tg = 0.6). More recently, the RS effect has been estimated, 
in the concordance model, for various sl ightly d i fferent 
valu es of the involved param eters (see, e.g. ICooravl ll2002h 
and IHu fc DodelsonI i2002h l: these estimations are not 
based on N-body simulations and they lead to sma ll values 
of Ar„ for £ ~ 700; thes e values are ~ 0.17 tiK in ICooravl 
Iggpj) and ~ 0.2 ^iK in IHu fc DodelsonI ll2002f) . As it has 
been discussed in section 4, the value of C700 extracted 
from our small RS simulations seems to be between 0.4 jiK 
and 0.5 fiK; this approximated value seems to be a little 
gr eater than those found (with othe r models and methods) 
bvlTuluie. Laeuna fc AnninosI l[l993) . ISeliakl (ll996l) . ICooravl 
J2002h . and lHu fc DodelsonI tOO'j) . The A^, spectrum given 
in the three last of these references and that obtained by 
us slowly decrease for £ > 700; nevert heless, o u r spec trum 
decreases slowe r than that shown in ICooravl ll2002ll and 
IHu fc DodelsonI (|20o3). Our A^ overestimations, at large £ 
values, associated to map smallness (see section 4) seem to 
explain this difference. 

The C3 correlation induced by t he R S ef- 
fect was studied bv iMollerach et alJ 1^9^ and 
iMunshi. Souradeep fc Starobingkil i^gsIT The main 
conclusion of these authors was that this odd correlation is 
smaller than its cosmic variance, which suggests that it is 
too small to be detected from a unique complete realization 
of the RS sky. Of course, we have not any map of the full 
sky, but very small maps of the RS effect and, consequently, 
it should not be possible the estimation of C3 from these 
small parts of the full sky. In spite of these considerations 
and taking into account that the mentioned estimates of C3 
are based on perturbation theory, which does not properly 
describe strong non-linear evolution inside cluster and 
substructures, we have tried the estimation of the C3 and 
C5 correlations from our 1.83° x 1.83° and 3.66° x 3.66° 
RS maps; nevertheless, this estimate has not been possible, 
which conflrms the low level of these odd correlations in the 
maps of the RS effect. 

Even correlations have been calculated for both CMB 
and RS small maps and, after comparisons, a certain de- 
viation with respect to Gaussianity has been pointed out 
in the RS case (see section 5); however, no deviations have 
appeared in the analysis of the superimposition of RS and 
CMB maps. It is due to the fact that the RS non-Gaussian 
part is much smaller than the CMB Gaussian one. The su- 
perimposition is almost Gaussian and, consequently, cosmic 
variance errors, ACi = [2/(2£ + 1)]^^^C^, are a good ap- 
proximation to the Ce uncertainty in full sky measurements. 
These errors can be easily calculated taking into account 
that the angular power spectrum of the superimposition 
is almost identical to that of pure CMB. In the £ interval 
[500,2000], where we have found our best estimate of the RS 
spectrum (see section 4), one easily finds the errors of A^, 
in jj,K, associated to cosmic variance. They are displayed in 
the solid line of Fig. 1131 In the same Figure, the dotted line 
is our best estimate of the RS spectrum (also represented 
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Figure 13. numbers in the vertical axis are temperatures in /iK. 
The solid Une displays the errors of due to cosmic variance 
and the dotted (dashed) line is the spectrum obtained from 
HR simulations for eg = 0.93 (erg = 0.76) 



in the dotted line of Fig. and the dashed Une is a new 
spectrum which has been obtained using the same method 
as in the case of the dotted line, but assuming erg = 0.76 
in our HR simulations. This new normalization is suggested 
by t he analysis of the three year WMAP data published 
iSpe rgcl ot al. 2006) after the submission of the first version 
of this paper. In Fig. 1131 we see that, for the two erg values 
we have considered, our upper limits to the RS spectrum 
are well below the cosmic errors for £ < 1200; hence, taking 
into account that WMAP cannot accurately measure Ce co- 
efficients for £ > 1200 (see Hin saw ct al.. (J006)), one easily 
concludes that the RS effect cannot be observed with this 
satellite. Finally, the upper limits to the RS spectrum dis- 
played in the dotted and dashed lines are well above 
errors for 1500 < £ < 2000; however, taking into account 
that the true RS spectra are expected to be around 0.2 jiK 
smaller than our upper limits for large £ values (see section 
4), it follows that the true values are not expected to 
be greater than the cosmic errors given by the solid line of 
Fig. 1131 This means that, although all the multipoles up 
to ^ ~ 2500 could be measured with PLANCK, the RS ef- 
fect would be too small to be detected with this satellite. 
All these considerations strongly suggest that the RS effect 
will not be directly detected in the near future. Indirect de- 
tection based on correlations of the CMB with the cluster 
distribution producing the effect deserves attention. 

Simulations in great boxes would lead to more extended 
RS maps (improving on power spectrum estimation and non- 
Gaussianity detection) . In the case of large simulation boxes 
with L ~ 1000 Mpc, which can be currently considered with 
appropriate N-body codes and computers, the magnification 
of the RS effect associated to periodicity is not negligible 
because eight of these boxes are yet necessary to cover the 
photon trajectory from z = 5.2 to present time. Of course, 
our method based on a preferred direction also applies for 
these big boxes. After minimum changes, our codes will be 
ready to built up better RS maps from appropriate N-body 
simulations, and also to statistically analyze these maps. 
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